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ABSTRACT 


The  present  paper  gives  an  algorithm  that  finds  simultaneously  the 
shortest  paths  between  many  pairs  of  nodes  in  a  given  network.  In  the  book 
by  Berge,  the  values  of  shortest  paths  between  many  pairs  of  nodes  are 
found.  Here,  we  use  a  special  matrix  multiplication  technique  to  find  the 
actual  arcs  that  are  used  to  form  the  shortest  paths.  In  a  network  with 
n  nodes,  log^fn-l]  special  matrix  multiplications  are  needed  to  find  all 
the  shorte-t  paths. 

The  present  paper  also  gives  an  algorithm  for  constructing  a  net- 

♦ 

work  with  prescribed  shortest  distances  and  with  the  total  distances  as¬ 
sociated  with  arcs  a  minimum. 


The  present  paper  gives  an  algorithm  that  finds  simultaneously  the 
shortest  paths  between  many  pairs  of  nodes  in  a  given  network,  an  algo¬ 
rithm  for  constructing  a  network  with  prescribed  shortest  path  distances, 
and  with  the  total  distances  associated  with  arcs  a  minimum. 


Many  papers  have  been  written  about  finding  one  shortest  path  from  a 
given  node  to  another  node  in  a  given  network  (see  for  example,  fl’j,  [3], 
[4],  [8],  or  the  review  [9]  ).  If  we  are  interested  in  shortest  paths  be¬ 
tween  many  pairs  of  nodes,  then  this  problem  can  certainly  be  solved  by 
using  any  of  the  existing  algorithm  for  finding  a  shortest  path  and  doing 
the  algorithm  over  and  over  for  every  pair  of  nodes.  Because  in  finding 
a  shortest  path  many  informations  are  obtained  during  the  process,  we  would 
expect  that  some  saving  would  be  acheived  if  we  want  to  find  many  shortest 
paths  Also,  for  finding  maximum  flow  values  between  many  pairs  of  nodes, 
it  is  known  that  a  great  amount  of  computation  can  be  saved  (see  [5]  ),  and 
the  known  dual  relationship  between  maximum  flow  and  shortest  path  in  an 
a-b  planar  network  would  also  lead  one  to  expect  that  many  shortest 
paths  can  be  found  simultaneously „  In  the  book  by  Berger  [2],  the  values 
of  shortest  paths  between  many  pairs  of  noeds  are  found  by  a  special  ma¬ 
trix  multiplication  technique;  in  Gomory  and  Hu  [0],  a  minimal  cost  net¬ 
work  i s  found  which  can  satisfy  given  flow  requirements  between  many  pairs 
of  nodes  simultaneously.. 

We  are  now  interested  in  finding  the  actual  arcs  that  are  used  to 
form  the  shortest  paths  Although  the  paper  is  self-contained,  we  shall 
use  many  of  the  notions  of  [2]  and  [6"„ 


We  consider  a  network  consisting  of  nodes  N  and  arcs  A  leading 

**■  J 
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from  nodes  to  ,  Associated  with  each  arc,  there  is  a  distance  func¬ 

tion  d  defined  on  it  We  shall  not  assume  either  d  «  d  or  that 

cl  c)  J  ^ 

diJ  +  ^Jk  ^  ^ik  *  but  we  a*5Ume  throughout  the  paper  that  distance 

functions  d  are  all  non- negative  . 

^  J 

We  shall  follow  many  authors  (see  for  example,  [2]  )  in  defining  two 
binary  operations 


(1) 


(2) 


Xi®k2  "  min  'X1  ’  X2* 
^l(*)^2  ”  ^1  +  ^2 


Using  this  special  addition  and  multiplication  instead  of  usual  ad¬ 
dition  and  multiplication,  we  have  a  special  matrix  product  of  the  two 
n  x  n  matrices  =  [d^]  anc*  ^ij"J  matrix  ^  3  [d^]  with 


J) 


^  »  minfd,  +  d'  }  s  =  1 

ij  is  sj 


n 


-2 


If  we  square  D  ^  [d  j  ,  the  distance  matrix  in  the  above  sense, 

-*•  J 


the  resulting  matrix  D1"  -*  [d"  j  would  give  the  values  of  shortest  paths 

iJ 

n-1 

using  two  or  less  arcs  between  any  two  nodes.  Similarly,  D  would 

give  the  shortest  path  values  between  any  two  nodes  using  n-1  or  less 

arcs  Since,  for  a  network  of  n  nodes,  shortest  paths  are  always  of 

,  %  PL 

n-1  arcs  or  less,  after  L  times  operations  xike  (3;,  D  Is  obtained 
with  2^^  n-1  and  L  <£  [log^n-l)  ]  where  [log^(n-l)  ]  indicates  the 
least  integer  greater  than  or  equal  to  log^n-l). 

In  the  paper  by  Gomory  and  Hu  76],  the  problem  is  to  find  a  mini¬ 
mum  cost  network  that  will  satisfy  all  given  flow  requirements  simul¬ 
taneously,  Assume  the  cost  of  building  an  arc  of  unit  capacity  between 

nodes  N,  ar.d  N  is  c  »  we  want  to  find  y  such  that 
i  J  i  J  1 J 
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J  i  p>(i 


(5)  #*SL  ^ 

i 


J  =*  P 


L  x”  *  yij 
p,q 


for  all  i  and  j 


(7)  min  z  =  cJ  .y.  .  where  R  are  given  flow  requirements  be- 
ij  iJ  PQ 

tween  nodes  and  .  This  problem  can  certainly  be  solved  by  using 
c  as  distance  and  then  find  the  shortest  paths  between  each  pair  of 

-L  J 

nodes  to  and  give  every  arc  in  the  shortest  path  the  amount  of 

y^j  =  R^^  .  The  final  y^  is  obtained  by  addition  ^  y^j  ■  y^  ^  .  In¬ 


stead  of  finding  shortest  paths  une  by  one,  we  can  find  the  cheapest  net¬ 
work  as  follows:  (see  [c]  ).  We  first  form  the  n  x  n  matrix  D  =»  [d  ] 

where  d  is  the  distance  from  nodes  N,  to  N  .  Then  we  form  the 
^  J  **■  J  , 

ok 

powers  of  D  by  squaring  in  the  sense  of  (3)  until  is  obtained 

where  2k  log0(n-l). 

2  4  8 

In  getting  D  ,  D  ,  D  ,  etc.,  we  need  only  keep  the  largest  power 

of  D  and  discard  the  previous  powers  of  D  .  At  the  same  time,  we  form 

matrices  B.  ,  B0  ,  where  B  »  [fc^J  and  b^ .  is  the  index  s  in 

]  ’  2  k  u  iJ  iJ 

(3)  for  which  the  minimum  was  obtained. 


..3. 


cl 

Aft‘-r  log^r.-l)  operations  or  less,  we  have  D'  and  ,  Bn  . 

Now  we  define  another  set  of  matrix  B  from  B^  .  We  start  with 

L,  -  RPq  ,  E  -  fb,  ,  1  is  obtained  by  adding  to  an  all-zero  matrix 
L+l  L  ij 

rP*  to  i,s  and  s,j  positions  where  ^  is  obtained  from  the  corresponding 

—  —  —  u+1 

Bt  .  Successive  B  are  defined  from  B  ,,  and  B  by  adding  b,  .  to 
L  u  u+1  u  ij 

an  all-zero  matrix  in  the  i;s,  and  s,j  positions  where  the  value  of  s  is 
obtained  from  B^  „  When  is  obtained,  gives  the  y  which  min 

(7)  and  that  will  satisfy  the  given  constraints  (4),  (5),  and  (6).  B.  ,  is 

Lt  _L 

identical  with  the  requirement  matrix  „  This  is  certainly  sufficient  to 

_url 

satisfy  all  requirements;  if  we  let  y^  =  b^  for  a11  *  and  J  • 

Since  we  want  to  find  the  cheapest  feasible  network,  we  try  to  find 
a  series  of  arcs  which  form  a  path  from  Ni  and  of  minimum  total 

cost.  B  is  obtained  by  adding  to  a  zero  matrix  to  the  two  pos- 

Ij 

itioas  in  B.  that  in  B  indicate  the  intermediate  node  in  which  the 

ij  Lt 

minimum  cost  is  achieved,  As  the  minimum  cost  path  for  each  R^  consists 


of  2L  arcs  or  less,  after  tracing  back  from  BT  .  L  times  to  B,  ,  B. 

LrM  1  1 


gives  the  y  which  form  the  cheapest  cost  feasible  network.  The  above 
algorithm  gives  the  network  which  contains  all  the  shortest  paths,  but 
does  not  give  any  information  abcut  whether  one  arc  is  used  in  a  given 
shortest  path  (or  requirement  RPq).  This  is  due  to  y  ^  being  the  sum 
of  many  requirements. 

To  convert,  the  above  algorithm  into  one  that  reveals  the  actual  arcs 
used  in  different  ,  we  need  only  artificially  give  rP^  as  follows. 


Let  R 
nodes 


,R.  ,  . , . ,R  be  the  given 
i  m  k 

We  shall  let  (8)  y  R^ 


m  requirements  between 
<  \.fl  for  k  ■=>  1, . . .  ,m-l 


m 


pairs  of 
For  example, 


i=»l 
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«  2k_1  will  satisfy  (8). 

—  k-1 

When  is  obtained,  we  can  uniquely  decompose  y  ■  Z6^2 

wher?  6^  «•  1  or  zero.  Suppose  we  want  to  find  the  shortest  path  from 

to  ,  first  find  ,  and  then  start  from  .  Following  the 

k-1 

arcs  with  2  in  the  decomposition  of  y^  will  lead  to  along 

the  shortest  path. 


Take  a  five-node  network  as  an  example.  The  distances  between  nodes 

2 

are  shown  in  Table  1.  (Note  distances  are  not  symmetric.  )  Hie  D  together 


with  matrix  is  shown  in  Table  2.  The  together  with  B ^  matrix  is  shown 
in  Table  3. 
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© 

© 

© 
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® 
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© 

B 

1 
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0 

5 

© 

B 

2 

2 

_ 

0 

© 

® 
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© 

0 

1 

3 

5 

5 

© 

B 

0 

2 

4 

B 

© 

B 

5 

0 

2 

2 

© 

2 

3 

5 

0 

0 

© 

2 

3 

5 

7 

0 

[di/  - 1)4 


Table  5 
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Suppose  that  we  want  to  find  the  shortest  paths 

N  to  N  ,  N  to  N  ,  N  ^  to  N  and  N  to 

1  4  2  5  3  b  ^ 

artificially  set  the  corresponding  flow  requirements 
8  ,  and  16  .  This  is  shown  in  Table  4. 


from  N  to  N 

1  } 

.  Then  we 


to  be  1,2,4 


© 

© 

© 

© 

© 

© 

0 

0 

1 

2 

0 

© 

0 

0 

0 

0 

4 

© 

0 

0 

0 

0 

8 

© 

16 

0 

0 

0 

0 

© 

0 

0 

0 

0 

0 

Table  4 


From  Bp  in  Table  3  and  in  Table  4,  we  get  Table  5  ( by  adding 


R54  to  i  ,  s  and  s  ,  j  )  . 


© 

© 

© 

© 

© 

© 

0 

2 

1 

0 

0 

© 

0 

0 

4 

2 

0 

© 

0 

0 

0 

0 

12 

© 

l6 

0 

0 

0 

0 

© 

_ 

0 

0 

0 

0 

_ 

0 

B 


r\ 

c. 


Table  p 
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From 


in  Table  2  and  in  Table  3,  we  aet  in  Table 


© 

© 

© 

© 

© 

0 

3 

0 

0 

0 

0 

0 

7 

0 

0 

© 

0 

0 

0 

14 

0 

© 

0 

0 

0 

0 

28 

© 

_ 

1  6 

0 

0 

0 

0 

B 


Table  6 


Note  that  every  number  in  Table  6  has  a  unique  decomposition  as 
Z  6  2  ,  For  example,  if  we  want  to  find  the  shortest  path  from 


k=»0 


fg 


Ng  to  N^  ,  then  we  look  through  and  find  out  what  numbers  contain 

4  as  their  partial  sum.  In  the  present  example,  we  have  7  a  4  +  2  +  1. 
14  a  8  +  4  +  2  and  l6  +  8  +  4  .  Therefore,  the  shortest  path  is 

from  N  to  N  ,  N  to  N,  ,  and  N,  to  N  .  If  numbers  are 

2  3  3  4  4  5 

represented  in  the  binary  system,  the  process  of  checking  whether 

k 

a  number  contains  2  is  a  very  easy  job. 

We  have  given  a  method  of  calculating  the  shortest  paths  between 
all  paire  of  nodes  in  a  network  in  which  distances  on  arcs  are  given. 

Now  we  consider  the  problems  of  realization  and  synthesis.  (The  problems 
of  realization  and  synthesis  of  undirected  networks  were  studied  by 
Hakimi  and  Yau. 
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Problem  of  realization:  Given  an  nxn  m  t  *  i>.  [£  ]  ,  Who. 

ij 

are  the  conditions  on  [£  ]  such  that  there  exists  an  n-node 

^  J 

network  whose  shortest  distances  arc  the  given  f  l.  .  ]  V  It  is  easily 

'  ij 

seen  that  the  necessary  and  sufficient  conditions  on  i  is  that 

J 

£  s  £  +  £  ,  for  all  i  ,  j  ,  k  .  The  necessity  is  clear.  If 

ik  ^  ij  jk 

£.,'>£,,  +  £.,  ,  then  the  shortest  path  from  N,  to  N  of 

ik  x  ij  jk  1  *  i  j 

distance  £  .  pi  is  the  shortest  path  from  N.  to  N  of  distance 
i.  *i  j  k 

£  is  a  path  from  I*.  to  N  and  of  less  distance  than  £  . 

jk  i  k  ik 

This  would  contradict  that  £  is  the  shortest  distance  from  N 
to  . 

The  sufficient  part  is  also  clear.  If  £J .  <  £J  <  +  £ ,,  for 

ik  iJ  jk 

all  i  ,  J  ,  k  ,  we  can  then  construct  n(n-l)  arcs  between  an 

n-node  network  with  d  .  =  £  for  all  i  ,  j  .  Then  the  shortest 

i.  J  X  J 

path  from  any  node  N  to  any  other  node  N.  is  just  the  single 

i  J 

fl¬ 
are  AJ  .  with  distance  d  .  .  Since  there  arc  many  n-node  networks 

ij  ij 

that  will  have  shortest  distances  equal  to  the  prescribed  ti  ]  > 

the  following  synthesis  problem  arises  naturally.  Given  an  nxn 

matrix  [£  ],find  a  n-node  network  with  given  [£  .]  as  shortest 
i-  j  i  1 

distances,  such  that  the  total  amount  of  d. .  associated  with  arcs 
'  ij 

of  the  network  is  minimum  (i.e.  min  S  d  )  .  From  now  on,  we  assume 

J 

dij  >  0  • 

We  shall  use  the  words  "optimum  network"  to  mean  the  network  with 
total  d  .  minimum  and  satisfying  the  prescribed  [£  ]  as  shortest 

J  J 

distances.  If  an  arc  A  is  in  an  optimum  network,  then  the  shortest 

ik 

path  from  N  to  N  consists  of  the  single  arc  A  ,  with 
1  I  k  ik 

d  -  £  .  For,  if  £  were  ?qual  tod  +d  +.,.+  d  <d  , 

ik  ik  '  ik  H  il  12  mk  ik  ’ 


the 


n  any  shortest  path  using  A  ,  can  use  the  arcs  A  A  ,...,  A 

iK  J  L  Id  J 


mk 


*  This  is  the  same  condition  given  by  Hakimi  and  Yau. 
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Instead  so  that  the  arc  A  can  b*  repla  ed  f  reduce  the  t  t  al 

IK 

Z  d  ^f  the  network.  This  weald  contradict  the  network  1°  an 

U 

optimum  network. 


Therefore,  If  an  arc  A  is  in  an  optimum  network,  then 

ik 

i  -i  d  ,  If  an  ar  A  is  net  in  an  optimum  network,  then 
ik  Ik  ik  ’ 


where 

from 


All ’  A12’* 
* 9 '  that 


are  arcs  of  the  optimum  network. 


It  foil  owe 


(10' 


£  >  max '  £  £  £ 

Ik  ^  il'  If;  '  mk 


From  '9),  the  ar:  ?  with  d  ,  e^ual  tc  the  smallest  £  must  be  in  an 

'  1J  i.j 

optimum  netv:rk„ 

If  w“  know  two  arcs  A,  A  a^e  in  any  optimum  network,  then 

ij  jk 

?  ~  £  -id  +  d  .  If  £  d  -t-  d  .  then  the  arc 

I*  S  Ij  jk  ij  jk  ik  ij  jk  ■ 

A  ,  will  not  be  in  any  optimum  network.  If  i  ,  <  d  f  d  , 
ik  ik  ij  jk 

then  A  mav  ''r  may  not  be  m  an  optimum  network.  This  depends 
ik 

cn  if  there  are  other  ar,_s  such  that  i  ,  "  d  r  d  .  However,  if 

:k  ip  pk 


d 

IJ 

<.  d 

cm 

and 

d  ,  <  d 
jk  ak. 

ior  any  nodes  m  ,  n 

d 

i.m 

+  d 

rrx 

for 

any  nede  m 

.  Consequently,  if 

•J 


f  d 


mk 


then  the  sr  -  A  !•;  in  an  optimum  network, 
ik 


From  the  above  onsideratiens,  we  :  an  develop  the  following 

ad.gon*hm  for  construe  ting  an  optimum  network.  Given  the  prescribed 

distance  matrix  [i  .J  ,  we  shall  circle  some  entries  £  and  pi* 

ij  ij 

bracket;  cn  other  entries. 
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Step  1 : 


Circle  the  smallest  £  .  which  has  not  yet  been  circled, 

bracketed,  or  crossed  cut.  v"  Intially,  none  of  £  ,  has  been 

J 

circled,  bracketed,  or  crossed  out.  )  When  an  £,  .  has  been 

ij 

circled,  cross  out  the  ith  column  and  the  jth  row.  Continue 

until  every  entry  is  circled  or  crossed  out. 

Step  2:  Form  the  matrix  D  with  d  .  equal  to  circled  or 

^  J 

bracketed  I,  ,  d  ,  -  ®  otherwise, 
ij  1J 

2  " 

Step  3:  Square  D  in  the  cense  of  (3).  Let  D  =  [d  ]  . 

tJ 

Compare  [d  with  [£  ]  .  If  £  is  not  circled  and  eaual 

ij  ij  ij 

to  d  "  ,  put  a  bracket  on  the  £j  .  . 

i  J 


If  every  £  .  is  either  circled  or  bracketed,  go  to 
i  j 


Step  4 . 

Step  3.  Otherwise,  erase  ail  crossed  out  columns  and  rows 
and  return  to  Step  1. 

Step  3- 

network 


Toe  matrix  with  d  •-  circled  £  .  is  the  optimum 

ij  iJ 


riote  in  Step  1  \£  1  is  being  circled  with  increasing  magnitude, 

ij' 

Because  of  the  cross-out  process,  the  [2  .  ]  circled  in  the 

1  vJ 

same  seep  is  not  connected.  Fr cm  (9)  and  (10)  they  must  be  in 

the  optimum  network.  In  step  3>  v-  heck  whether  or  not  some  of 

these  circled  arcs  form  a  path  and  if  they  would  imply  the  orussion 

of  some  other  ar^s.  Since  we  only  square  the  distance  matrix  formed 

in  Step  2,  we  only  check  paths  coasistlng  of  two  arcs.  When  an 

entry  (£  I  is  equal  to  d  +  d  ,  this  £  is  bracketed 
"  ik  ij  Jk  lk 

and  will  be  used  tc  form  the  matrix  D  in  Step  2  of  the  next  cycle. 

As  at  the  end  cf  Step  1  we  erase  fill  crossed  out  columns  and  rows, 

every  entry  £  will  sooner  •- r  later  by  circled  -or  bracketed  which 
1  k 


implies  that  arc  is  in  or  is  not  in  the  optimum  network.  The 
uniqueness  of  the  optimum  network  follows  from  the  observation  that 
there  is  no  alternate  choice  of  circling  an  arc. 
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